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95: 1. Introduction 



The peculiar properties of the unstable particles (UP) and resonances were being discussed 
during the last decades. Among them, the assumption that the decay of UP or resonance 
(R) proceeds independently of its production remains of interest [H El E]. Formally, this 
effect is expressed as the factorization of a cross-section or decay rate [3]. The processes of 
type ab — » Rx cdx were considered in Ref. [3j. It was shown, that the factorization always 
is valid for a scalar R and does not take place for a vector and spinor R. The factorization 
usually is related with the narrow-width approximation (NWA) [4J, which makes five critical 
assumptions [5]. 

There is another way to get factorization effect, which is connected with propagator 
structure [6J. The decay processes of type a — > Rx cdx, where R is UP with a large width, 
were analyzed systematically in Ref. [6]. It was shown in this work, that the factorization 
always is valid for a scalar i?, while for a vector and spinor R it occurs when the propagators' 
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numerators are ri^^{q) = g^^ — q^qv/q^ and 57(g) = q + q, respectively, where q = qiY and 
Q = V QiQ^- Such a structure of propagators always provides the exact factorization for any 
tree process and is an analog of NWA, which is discussed in Section 3. These propagators 
were constructed in the model of UP with a smeared mass [7J and describe some effective 
(dressed by self-energy insertion) unstable fields. Note that the structure of the expressions 
77^1/ (g) and 17(5) is not related with the choice of the gauge (see the third section). 

In this work, we systematically analyze the processes of type ab R cd, where R is 
scalar, vector or spinor UP with a large width (or resonance) and a, b, c, d are the stable or 
long-lived particles of any kind. It was shown that the cross-section a{ab — > i? ^ cd) can be 
represented in the universal factorized form when the same expressions ri^y{q) = g^v — QuQv/ (f 
and fi{q) = q + q are used to describe the propagator's numerator of vector and spinor UP, 
respectively. This result have been received strictly by direct calculations for all types of 
particles a,b,c,d and R (Section 2). The factorization approach is applied in Section 3 for 
the complicate processes of scattering with the consequent decays of the final states. In 
Section 4, we analyze some methodological and phenomenological aspects of factorization. 

2. Universal factorized formula for the 
cross-section of two-particle scattering 

In this section, we consider inelastic scattering of type a6 — > i? — > cd, where R is the UP 
with a large width in s-channel and a, b, c, d are stable (quasi-stable) particles of any kind. 
The vertexes are defined by the Lagrangian in the simplest standard form: 

^fc =5-00102; 5'0^iV'2; g(j)Vif,V2^; 5(V^(0f02 - 02'0i); gVf^tPiYicv + CAi5)i'2; 

gV,,V2.V^[g^''{p2-PiT + g'''{2pi+P2r - g'^'ipi + 2^2)^']. (1) 

In the expressions (II]) 0, and ip are the scalar, vector and spinor fields, respectively, pi and 
P2 are the momenta of the particles a and b (or c and d). 

Here we show, that the cross-section a{ab cd) can be expressed in a factorized 

universal form in terms of decay widths r(i? — > ab) and r(i? — > cd), when the expressions for 
propagators' numerators rj^y{q) = g^^, — q^qu/q^ and 17(g) = g + g are used. This expressions 
are constructed within the model of unstable particles with a smeared mass, which briefly 
considered in Appendix. The vahdity of these expressions have been discussed in Refs. [HI 
|7] and will be considered in the third section. It is convenient to employ the universal 
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expressions for widths r(_R — > ab) and T{R cd) in a stable particle approximation [6]: 

Ti{R a6) = — A(ma, m,,; mR)fi{ma, mb] mRj, (2) 



where m|. = g^, = (pi + ^2)^ and: 

(3) 



A(m,, m,; m^) = [1 - 2 " , ^ + ^^^_^]i/2. 



2 4 

The same expressions and relations are in order for the width r(i? — > cd). The functions 
fiinia, mb] mn) are defined by the corresponding vertexes. If these vertexes are described by 
Eqs.([T]), then the functions (further we omit the arguments) in tree approximation are 
defined by the following expressions [6J: 

1 1 (m'i — m? — ml)'^ 

^ - M2, h = ^; - ^1^2, /2 = —[1 + \ I 2 ]; 

2m0 Smfm2 
^^0102, /5 = ^A^(mi,m2;my); V^Vicf), /g = -^[1+ 

O oTfly 

{m\r + m\- 

~^ o 2 2 J ' 

8mymf 

- 2 mf + ml (m2 + m2)2 mims 
^ = 3"^^^^+[' - ^< 2<^] + 



1^ ^ ^iV^2, /s = 7T7ri^[l + + /W2) - 2(9/i? + 16/ii/i2 + + 8(/i?- 



rrr( 
24mf m2 

.2, .2 



4/ii/i2 - 4/ii/i2 + /i2) + /^i + 8aH/^2 - 18Aii/i2 + 8/ii/i2 + /i2], /^i,2 = 2My; 

^ 0^^!, /g = ^1 + 2 — + — 2 — ; 
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^ m^-m^)^ ^ m^ + m^ -2m^ - 3c ^ 
2m2 2m^ ' ^ 



^ ^ ^V'l, /lo = m^{c+[ ^^2^2 + ^::;:2 ] - 3c-— }; 

2 I 2 2 2 / A\ 

c+ = Cy + c^, c_ = Cy - . (4) 



Note that the function /§, given in Ref. [6], contains an error and we give here corrected ex- 
pression for this function. It is convenient in the further calculations to employ the relations, 
which take place in the center-of-mass system: 

1 22 1 22 

n 1 r — fnt T n 1 r fni — m„ , 

p\ = 2?[1 + ^^]' P2 = 2^[1 + 

(PiQ) = \{<f + ^1- ^1), {P2q) = ^{q^ + ml), 

(Pm) = ^{q'^ -ml-ml), \pi\ = \p2\ = ^qXma, mb, q) . (5) 



3 



The analogous relations occur for the momenta ki and k2 of the particles c and d. In Eqs.(l5]) 
the symbol q has different meanings in the expressions (piq), Q = Pi + P2 (q is 4-momentum) 
and in the expression g[l + f{q)], where q = \J{q- q) is a number. 

With the help of the relations (I2l)-(l5]) and above discussed expressions for propagators, 
we have got by tedious but straightforward calculations the universal factorized cross-section 
for all permissible combinations of particles (a, 6, R, c, d) : 

a(ab^R^cd) = lQ-i2Jn + l) rg(.)rg(.) 

In Eq.([6]) Jfc is spin of the particle {k = a,b, R), s = {pi + P2Y , r^''(s) = V{R{s) ah) and 
Pr{s) is propagator's denominator of the UP or resonance R. The expressions for r^^(s) and 
r^(s) follow from Eqs.(l2l|4l), when squared mass of UP is m\ = q^ = s. The factorization 
of cross-section does not depend on the definition of Pr{s), which can be determined in 
a phenomenological way, in Breit-Wigner or pole form [H [9], etc. The expression (|6]) is 
a natural generalization of the spin-averaged Breit-Wigner (non-relativistic) cross-section, 
defined by the expression (37.51) in Ref. |10j. Note that the factorization is exact in our 
approach, while in the traditional one it occurs as an approximation. 

The cross-section of exclusive process ah R cd, defined by Eq.Q, does not depend 
on Jc and Jd- So, it can be summarized over final channels R — > cd: 

In Eq.([7]) ki = 2Jj + 1 and r^*(s) = ^c(^r^(s), where for simplicity we restrict ourselves by 
two-particle channels. 

The factorization effect, expressed by Eq.(l6]), has two aspects. On the one hand it means 
that the decay of UP proceeds independently of its production in the approach considered. 
On the other hand it leads to significant simplification of calculations, in particular, in the 
case of the complicate processes (the scattering with chain decay of products). 



3. Cross-section of the process ah ^ R ^ RiX cdx 

Here, we consider factorization effects in the case of complicate chain processes. For example, 
let us discuss the process of scattering ah ^ R ^ Rix with consequent decay Ri cd. In 
this case, Eq. ([6]) has the form: 
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To calculate the value r^^^(s) we apply convolution formula, which accounts FWE in the 
decay R{s) Rix |6]: 



r{R{s)^Rix)= / T{R{s)^ Riiq)x)pR,{q)dq. (9) 

In Eq. q = PR—Px, <ii,2 are defined by kinematics of the process and pR^q) = gr^*(g)/7r|PR(g)p 
is interpreted in the model [6] as distribution function of the smeared mass of unstable 
particle Ri. Convolution structure of Eq. (Q is caused by factorization of decay rate 
r(_R — i> Rix x,aU). This effect takes place exactly when the model propagators 
and rj^y{q) are used (as in the present work). 
From Eqs. ^ and ^ it follows: 



a{ah ^ R ^ Rix) 

^2 m{s)^Ri{q)x)pRAq)dq'. (10) 



kahX^{ma,mi,]s/s)\PR 
Using the expression for PR^{q), from Eq. flTUl) we can get the cross-section of exclusive 
process, for example ab R Rix — > cdx. To this effect we represent r^*(g) in the form: 

^rM) = J2^rM ^rM = nRiiq) - cd) . (11) 

As a result, from ( ITOi) and ( fTTl) we get: 



a{ab —>■ R ^ cdx) = 

_ IQkR rg(g ) ^^V.o..^ , o ^rff^(^) 

kahX^{ma,mi,;y/s) |Pr,(s) 
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r(i?(.) ^ i?i(g)x)^p|^ rfg^ (12) 



Similar structure arises in the case R — > R1R2, i.e. when there are two UP in the 
final state, which have two-particle decay channels (semi-analytical approach). Thus, the 
model gives a convenient instrument to describe two-particle scattering accompanied by 
complicated decay-chain processes. However, we have checked by direct calculations only 
two types of processes - the decay of type a Rx hex [6\ and the scattering of type 
ab R cd. The more complicated processes, such as decay a —>■ R1R2 —>■ cdef and 
scattering ab ^ R ^ R1R2 cdef, will be the subject of the next paper. 

4. Methodological and phenomenological 
analysis of the factorization effect 



The model factorization of a decay width and cross-section of the processes with UP in an 
intermediate state was established by straightforward calculations at tree level. However, 
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these calculations in the effective theory of UP [7] account for some loop diagrams. The ver- 
tex and self-energy type corrections can be included into r^(s) and Pr{s) respectively. These 
corrections do not breakdown a factorization, but the interaction between initial and final 
states does. However, such an interaction has no clear and explicit status in perturbation 
theory due to UP (or resonance) is not a perturbative object in the resonance neighbor- 
hood [71 [m [12]. As it was noted in Ref. |l3j, such non-factorizable corrections give small 
contribution to the processes e~^e^ ZZ, WW, 4/ near the resonance range. 

Now we consider another aspect of factorization effect, namely, the determination of 
dressed propagator of UP. Factorization of decay width and cross-section does not depend 
on the structure of propagator's denominator PR^q), but crucially depends on the structure of 
its numerator in the case of vector and spinor UP. As it was verified by direct calculations, 
the factorization always takes place in the case of scalar UP. The expressions Ti^y{mR) = 
OfMu — QfiQu/^'n and fj{mji) = q + tjir for vector and spinor UP, respectively, do not lead to 
exact factorization. But the expressions rj^y{q) = g^^ — <i^l(lu/(f and 17(g) = q + q strictly lead 
to factorization for any kinds of other particles. It should be noted that the definition of the 
functions rj^y{q) and fj{q) is not related with the choice of the gauge, because effective theory 
of UP [7] is not the gauge theory. The choice of q instead of mn in the rj^^j and 17 may seems 
contradict to the equation of motion for vector and spinor UP. However, this statement is 
valid for the stable particle with fixed mass. In the case of UP the question arises what the 
mass participates in equation of motion - pole mass or one of the renormalized mass [14j? 
An account of uncertainty relation by smearing of mass intensifies the question. There is 
no unique and strict determination of dressed propagator structure for vector and spinor 
UP due to the specific nature of Dyson summation in these cases ^ . The situation is more 
complicated and involved in the case of hadron resonance. So, the functions rj^^, and 17 have 
rather phenomenological (or model) than theoretical status. The model of UP [7] defines 
these functions as ri^^{q) and 57(g), which describe the dressed propagators of UP in the 
resonance neighborhood. 

Further, we briefly analyze the phenomenological aspect of factorization. In the low- 
energy experiments of type e~^e~ p,uj... — * tt+tt^,... we can not distinguish propagators 
Vfj-vi^R) and rj^^{q) even for the wide resonance. This is due to the equality e~(pi)(]3i + 
V2)^^{p2) = 0, when the functions rj^i, reduce to g^^ in both cases. In the high-energy 
experiments of type e"'"e^ — > Z ^ //, where / is quark or lepton (we neglect 7 — Z 
interference), the transverse part of amplitude is 

M, ~ e^(pi)g(ce - 75)e-(p2)/+(A;i)g(c; - (13) 
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where q = P1+P2 = ki + k2. From Eg. (1131) with the help of the Dirac equations in momentum 
representation it follows 

Mg ~ mem/e-(pi)75e-(p2)/+(fci)75/+(fc2). (14) 

As a result, we get the terms meimf/q^ and meimf /m?^ for 77(g) and ri{mi{), respectively. The 
difference of these values is of the order of {mefrtf /m^^) ' i^z — (l)/^z at energy ~ m|. 
Thus, the distinction between the structure of two type of the expressions r/^j/ is negligible 
in a wide range of energy. 

The structure of 17 can be studied in the process of type VF R ^ V' F' , where V and 
F are vector and fermion field, R is, for instance, baryon resonance with a large width. In 
this case, the difference between fj{mpi) and 57(g) is characterized by the value ~ Tn/mFt at 
peak region, and this problem demands more detailed analysis. 

From this analysis it follows that method of factorization is a simple analytical analog of 
narrow- width approximation (NWA, which contains five critical assumptions [5]). Instead, 
we use the structure of propagators' numerators //(g), which follows from usual ones under 
a simple transformation m/j g, and one assumption: there is no significant interference 
with non- resonant processes (fifth assumption of NWA). The rest assumptions of NWA can 
be derived from the first our point, where some of them are not obligatory in the special 
cases. The method leads to factorization in two type of processes - in the decay-chain [6] 
(universal convolution formula) and scattering ones (universal formula ([6])). Combining these 
two results, we get a simple and strict algorithm of analytical description of the complicated 
processes. 

5. Conclusion 

The factorization effect gives us a convenient phenomenological way to describe the three- 
particle decays and two-particle scattering processes. This effect significantly simplifies cal- 
culations and gives compact universal formulae for the decay rate and cross-section. 

In this work, we have shown that the factorization always is valid when scalar UP is in 
the intermediate state. In the case of vector or spinor intermediate states, the factorization 
takes place when the specific propagators are used for these states. These propagators are 
derived in the model of UP with a random (smeared) mass. They negligibly differ from the 
traditional propagators at peak area and follow from the smearing of mass in accordance with 
the uncertainty relation. Our method makes it possible significantly simplify the calculation 
of the complicated decay-chain and scattering processes. It is some analytical analog of NWA 
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and gives a simple and strict algorithm for calculations. This approach can be treated also as 
convenient approximation, which always is valid in the resonance range, where non-resonance 
contribution is small. 

We have fulfilled also a short methodological and phenomenological analysis of the ap- 
proach under discussion. It was shown, that in the process e~^e~ // the difference between 
two forms of propagators is negligible in a wide range of energy. It can be significant in the 
processes with baryon resonance in an intermediate state, but in this case we should fulfill 
an additional analysis. 



6. Appendix 

In this section, we briefly describe the model of UP with a smeared mass and construct the 
propagators for the vector and spinor fields. The structure of these propagators lead to the 
factorization effect in the processes with the participation of the UP in the intermediate 
state. The model field wave function, which describes UP, is represented in the form [7]: 

^aix) = j (15) 

where $a(x,yu) is standard spectral component, which defines a particle with a fixed mass 
squared = /i in the stable particle approximation (SPA). The weight function is 
formed by the self-energy interactions of UP with vacuum fluctuations and decay products. 
This function describes the smeared (fuzzed) mass-shell of UP. 

The model Lagrangian, which determines a "free" (effective) unstable field $(x), has the 
convolution form: 

L($(x)) = j L{^{x,^^))\^{^^)\^d^^. (16) 

In Eg. (1161) L($(x,/i)) is the standard Lagrangian, which describes model "free" field com- 
ponent $(x, /i) in the stable particle approximation {m? = fi). 

From Eq.( fT6|) and prescription c}$(x, /i)/(9$(x, fi') = 6{fi — fi) it follows the Klein-Gordon 
equation for the spectral component of scalar or vector field: 

(□-/i)$,(x,/i) =0. (17) 

In analogy with f|T7|) one can get the Dirac equation for fermion spectral component. As a 
result, we get the standard representation of the field function $a(x, /i) with a fixed mass 
parameter /x (spectral component). All standard definitions, relations and frequency expan- 
sion take place for $a(A;, /x), however, the relation /c° = \/k^~+~fl defines the smeared (fuzzy) 
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mass-shell due to a random nature of the mass parameter /i. The convolution (diagonal) 
representation of the "free" Lagrangian f|T6l) has an assumption (or approximation?) that 
the states with different fi do not interact in the approximation of the model "free" fields. 

The expressions f ll5l) -(|T7 l) define the model "free" unstable field as some effective field. 
As it was mentioned above, this field is formed by an interaction of "bare" UP with the 
vacuum fluctuations and decay products, that is includes self-energy contribution in the 
resonant region. Such an interaction leads to the spreading (smearing) of mass, described by 
the function uj{fi) or p{fi) = \uj{fi)\'^. Thus, we go from the distribution p^^{fi) = S{fi — M^) 
for "bare" particles to some smooth density function p(/i) = |ci;(/i)p with mean value fl ~ 
and mean square deviation o"^ ^ F. So, the UP is characterized by the weight function uj{fi) 
or probability density p{p) with parameters M and F (or real and imaginary parts of a pole). 

The commutative relations for the model operators have an additional 5-function: 

[Kik,p), $^(g,/i')]± = 6{p- fi')6{k-q)6o,p, (18) 

where subscripts ± correspond to the fermion and boson fields. The presence oi 5{p — p) 
in Eq.f ll8p means an assumption - the acts of creation and annihilation of the particles with 
various /i (the random mass squared) do not interfere. Thus, the parameter p has the status 
of physically distinguishable value of a random m^. This assumption is naturally related 
with a diagonal form of Eqs.f ll6l) and (ITTl) and directly follows from the interpretation of 
as a random parameter /i. By integrating the both sides of Eq. ffTSj) with weights uj*{p)uj{p ) 
one can get the standard commutative relations 

[K{k)Mpm± = ^{k-q)5^p, (19) 

where (A;) is the full operator field function in the momentum representation 

^a{k) = j ^^ikf^Mf^)d^i. (20) 

It should be noted that Eq.(l20l) follows from Eq.( |T9i) when f \u{p)\'^dp = 1, that is |c<j(yu)p 
can be interpreted as a normalized probability density. 

The expressions ( IT5|) . (fT6l) and (fTS!) are the principal elements of the model. The weight 
function uj{p) (or p{p)) is full characteristic of UP in the framework of the model. The 
relations (|T8l) define the structure of the model amplitude and transition probability. 

Here, we consider the model amplitude for the simplest processes with UP in an initial 
or final state and get the convolution formula as a direct consequence of the model. The 
expression for a scalar operator field [7] is 

= / J ^tII^^'"^^^"' (21) 
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where g° = + /i and a^{q,^) are the creation or annihilation operators of UP with the 
momentum q and mass squared rn? = /x. Taking into account Eq. (|T8|) one can get: 

[d-(fc,/i),0+(x)]_; r(x),a+(fc,/i)]-= .^^^^ e^''^ (22) 



where = \/k^~^-~fi. The expressions fl22p differ from the standard ones by the factor uj{fi) 
only. From this result it follows that, if a~^{k, fi)\0) and {0\d~{k, fi) define UP with the mass 
m = y/JI and momentum k in the initial or final states, then the amplitude for the transition 
$ 001 is 

A{k,iJ,) =uj{fi)A'\k,ix), (23) 

where A^^{k,fi) is the amplitude in the stable particle approximation. This amplitude is 
calculated in the standard way and can include the higher corrections. Moreover, it can be 
an effective amplitude for the processes with hadron participation. From Eq. fl23p it follows 
that the differential (on /i) probability of transition is dP{k,fi) = p{fi)\A{k, ^)\'^dfi. 

To define the transition probability of the process $ 00i, where is UP with a 
large width, we should take into account the status of the parameter /i as a physically 
distinguishable value, which follows from Eq.(|T8l). Thus, the differential (on k) probability 
is 

dT{k) = j dr\k,fx)p{fx)dfx. (24) 

In Eq.(!24l) the differential probability dV^^^k^ii) is defined in the standard way (the stable 
particle approximation): 

dr\k, p) = ^S{k^ -ks- k^)\A'\k, ii)\^dksdh , (25) 

where k = (/c$, k^j,, ki) denotes the momenta of particles. From Eqs.(l2^ and (l25l) it directly 
follows the well-known convolution formula for a decay rate 

r(m$,mi)= / r''*(m$, mi; /i)p(/i)(i/i , (26) 
J /ii 

where p(/i) = |ci;(/i)p, fii and fi2 are the threshold and maximal invariant mass squared of 
an unstable particle 0. 

An account of higher corrections in the amplitude fl2Sl) keeps the convolution form of 
Eq. fl26l) . This form can be destroyed by the interaction between the products of UP (0) decay 
and initial $ or final 0i states. The calculation in this case can be performed in the standard 
way, but UP in the intermediate state is described by the model propagator. However, the 
calculation within the framework of perturbative theory (PT) can not be applicable to the 
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UP with a large width, that is to the short-hving particle. In any case, the applicability of 
the PT, of the model approach or convolution method to the decays considered should be 
justified by an experiment. The validity of the CM was demonstrated for many processes, 
but this problem needs in more detailed investigation. If there are two UP with large widths 
in a final state $ 4'i4'2, then in analogy with the previous case one can get the double 
convolution formula: 

T{m^) = j j T'\m^; fii,fi2)Pi{fJ'i)P2{fJ'2)dfiidfi2- (27) 

The derivation of CF for the cases when there is a vector or spinor UP in the final state 
can be done in analogy with the case of scalar UP. However, in Eqs.l pTl) . ( l22l) and ( !23l) one 
should take into account the polarization vector em{q) or spinor u'^'^{q), where momentum 
q is on fuzzy mass-shell. As a result, we get the polarization matrix with = /i. In the 
case of vector UP in the final state we have 

^rn{q)el{q) = -gmn + grnQn/P , (28) 

e 

and in the case of spinor UP in the final state: 

^<'±(g)«;'^(g) = ^(q^^)^^, (29) 

where the summation over polarization is implied and g° = ^/q^ + p- The same relations 
take place for the initial states, however one have to average over the polarizations. The 
formulae fl26l) and fl27|) describe FWE in full analogy with the phenomenological convolution 
method. Similar method, called the semi-analytical approach, was applied in calculations 
of cross-section of the processes e~^e^ ZZ,WW at LEP2 energy where the phase 
space of the final states was integrated in analogy with fl27j) . This approach gives a simple 
expressions for the cross-sections, which are equivalent to the inclusive cross-section of the 
fourth-fermion reactions. The calculation of this cross-section in the framework of standard 
PT is very complicated and usually carried out with the help of the Monte-Carlo simulation. 
The model under consideration gives a quantum field basis for CM, which takes into account 
the fundamental uncertainty relation, provides a simple expressions for decay rates and is in 
a good agreement with the experimental data on some processes. To evaluate FWE in the 
case, when UP is in an initial state, we have to take into consideration the process of UP 
production. If UP is in an intermediate state, then the description of FWE is equivalent to 
the traditional one, but the propagators are determined by the model. 

Now, we consider the the structure of the model propagators. With the help of the 
traditional method, one can get from Eqs. ffTSl) . f|T8|) and (|20|) the expression for the unstable 
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scalar Green function [7]: 

{Q\T{ct>{x),cl>{y))\Q) = D{x-y)= I D{x ~ y, fi)p{fi)dfi . (30) 



In Eg. (1301) D{x,fi) is a standard scalar Green function with = /i, which describes UP in 
an intermediate state: 



D{x, fi) = -— / ———--dk . (31) 
(zvr)* J — fi + te 

The right-hand side of Eq. fl30|) is the Lehmann-like spectral (on fi) representation of the 
scalar Green function. Taking into account the relation between scalar and vector Green 
functions, we can get the Green function of the vector unstable field in the form: 

1 

9^n-k^kn/^l^_ikx^^_ (32) 



(27r)4 J k^ - j_M + ie 
Analogously, the Green function of the spinor unstable field is 

D{x, /.) = {^^ + ^)D{x, /,) = ^ I -^±^^e-'^^dk , (33) 

where k = kiY- These Green functions in momentum representation have a convolution 
form: 

Dmnik) = J Dmnik,n)p{n)dn, D{k) = J D{k, fi) p{ij,)dn . (34) 

The expression fl32tl34p for the propagators of vector and spinor fields leads to the effect of 
factorization, which in turn, gives the convolution formula. 
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